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ABSTRACT 

The National Combustion Code (NCC) was 
used to calculate the steady state, non-reacting flow 
field of a prototype Lean Direct Injection (LDI) swirler. 
This configuration used nine groups of eight holes 
drilled at a thirty-five degree angle to induce swirl. 
These nine groups created swirl in the same direction, 
or a co rotating pattern. The static pressure drop across 
the holes was fixed at approximately four percent. 
Computations were performed on one quarter of the 
geometry, because the geometry is considered 
rotationally periodic every ninety degrees. The final 
computational grid used was approximately 2.26 
million tetrahedral cells, and a cubic non-linear k - 
epsilon model was used to model turbulence. The NCC 
results were then compared to time averaged Laser 
Doppler Velocimetry (LDV) data. The LDV 
measurements were performed on the full geometry, but 
four ninths of the geometry was measured. One, two, 
and three-dimensional representations of both flow 
fields are presented. The NCC computations compare 
both qualitatively and quantitatively well to the LDV 
data, but differences exist downstream. The comparison 
is encouraging, and shows that NCC can be used for 
future injector design studies. To improve the flow 
prediction accuracy of turbulent, three-dimensional, 
recirculating flow fields with the NCC, 
recommendations are given. 

INTRODUCTION 

Most combustion processes have, in some way 
or another, a recirculating flow field. This recirculation 
stabilizes the reaction zone, or flame, but an 
unnecessarily large recirculation zone can result in high 
NOx values for combustion systems. The size of this 
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recirculation zone is crucial to the performance of state 
of the art, low emissions hardware. If this is a large- 
scale combustion process, the flow field will probably 
be turbulent, and therefore three-dimensional. This 
paper will deal primarily with flow fields resulting from 
Lean Direct Injection (LDI) concepts. LDI is a concept 
that depends heavily on the design of the swirler. The 
LDI concept has the potential to reduce NOx values 
from 50 - 70 % of current values, with good flame 
stability characteristics. It is cost effective and 
(hopefully) beneficial to do most of the design work for 
an LDI swirler using Computer Aided Design (CAD) 
and Computer Aided Engineering (CAE) tools. 
Computational Fluid Dynamics (CFD) codes are a CAE 
tool that can calculate three-dimensional flows in 
complex geometries. However, CFD codes are only 
beginning to correctly calculate the flow fields for 
complex devices, and the related combustion models 
usually remove a large portion of the flow physics. 

The National Combustion Code (NCC) is a 
state of the art CFD program specifically designed for 
combustion processes. A short summary of the features 
of NCC pertaining to this paper are: the use of 
unstructured grids, 1 massively parallel computing - 
with almost perfectly linear scalability, 2 a dynamic wall 
function with the effect of adverse pressure gradient, 3 
low Reynolds number wall treatment, 4 and a cubic non- 
linear k-epsilon turbulence model. 5, 6 The combination 
of these features is usually not available in other CFD 
codes and gives the NCC an advantage when 
computing recirculating turbulent flows. These features 
need to be validated, before the NCC is accepted as a 
design tool. The NCC has previously undergone 
benchmarking for simple flows, 7 and large-scale 
validations are being conducted. 8 

The purpose of this study is to quantify how 
well the NCC calculates a turbulent, three-dimensional, 
recirculating flow field. The comparison is against 
three-dimensional Faser Doppler Velocimetry (FDV) 
measurements on a prototype EDI swirler. This 
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comparison should show that the NCC could be a 
useful design tool for LDI injectors. 

THE LDI BENCMARK CASE 

Geometry and Computational Grid 

The LDI swirler is comprised of holes drilled 
at a 35 degree angle to induce swirl. There are nine sets 
of eight holes rotating in the same direction, or co 
rotating. Figure 1 and Figure 2 show the LDI swirler 
solid geometry at various angles. The swirler was 
placed in a 76.2 mm square pipe. (As indicated by a 
square box placed around the swirler.) The solid 
geometry model was created using Pro Engineer. The 
computational grid was generated using the CFD- 
GEOM program. 9 A tetrahedral grid was used 
throughout the computational domain. Hybrid 
prismatic/tetrahedral grids were attempted, but the 
quality of the prism cells was unacceptable. Looking at 
Figure 1 and Figure 2, it is obvious (to the author) that 
the geometry is rotationally periodic every ninety 
degrees. A special modification was made to the model 
because of the periodic boundary conditions. A small 
“sting” was added to the centerline of the Pro Engineer 
model. This “sting” was added to prevent both periodic 
boundary conditions from touching one cell. This 
condition is not desirable, because it causes numerical 
instabilities. (This is because the NCC periodic 
boundary conditions can have a different number of 
cells on each side.) To create this section, the Pro 
Engineer model was modified using the “cut” 
operation. The Pro Engineer model was then translated 
into IGES format, for use as Non Uniform Rational B - 
Splines (NURBS) surfaces in CFD-GEOM. This 
method was very time efficient. It took only 20 
minutes to generate the largest meshes using this grid 
generation “system”. During the start of this study, 
computational grids with 60,000 cells were used. The 
grid was refined until acceptable results were obtained. 
The final computational grid used was approximately 
2.26 million cells. Figure 3 and Figure 4 show the 
computational grid for the entire domain and at the 
swirler. 



Figure 1 - Rear View of the LDI 35 Degree Co 
Rotating Swirler 



Figure 2 - Front View of the LDI 35 Degree Co 
Rotating Swirler 



Figure 3 - Overview of the Final NCC 
Computational Grid at 2.26 million cells 
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L = 190mm 



Figure 4 - Close up of the final NCC 
Computational Grid at the Swirler Region 


Flow Conditions 

The LDI 35 degree co rotating swirler 
geometry was only computed as a non-reacting flow 
case. The NCC computations and LDV measurements 
were performed at approximately atmospheric pressure, 
with a 4 % static pressure drop across the swirler. 
Figure 5 shows an overview of the flow conditions for 
both the NCC computations and LDV measurements. 
The turbulence intensity was assumed at 10 % and the 
length scale was assumed 3% the inlet diameter. The 
exit pressure was atmospheric, and the exit temperature 
was at 298 K. The flow was confined with windowed 
walls for optical access. 

NCC Computations 

The NCC computations were run in general 
until the residuals were reduced three orders of 
magnitude. The mass flow rates at the boundary 
conditions were also monitored as a convergence 
criterion. Dissipation was set at .05 for second order 
dissipation (e2) and .08 for fourth order dissipation 
(84). 10 After the residual dropped three orders of 
magnitude, the 82 was set to zero, and the s4 was set to 
0.03125. The value of k2, the constant that scales the 
second order dissipation pressure switch, was set at 
0.25. Setting the second order dissipation to zero 
greatly improves the numerical accuracy, but if this 
initialization procedure was not done, the NCC 
sometimes ran into stability problems. A CFL number 
of 2.0 was used. 

The cubic non-linear k-epsilon model was 
used. This model was chosen because of its superior 
performance for computing swirling flows, while being 
computationally efficient compared to Reynolds stress 
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Figure 5 - Overview of Boundary Conditions for 

the LDI 35 Degree Co Rotating Swirler 

turbulence models. No other turbulence models were 
used or compared. Various turbulent wall conditions 
were used during the investigation. A “standard” wall 
function was originally used, but this gave extremely 
poor results. This wall function was then replaced in 
the NCC with a dynamic wall function with pressure 
gradient effects. This wall condition, along with the 
low Reynolds number variant of the cubic non-linear k- 
epsilon model, gave very similar results. Because this 
wall function is so new, and not well validated, the low 
Reynolds number wall treatment results are reported. 

Computations were originally performed on a 
variety of computer platforms, namely SGI Origin 
2000’s at NASA Ames and Linux clusters at NASA 
Glenn. Only Origin 2000 ’s performed calculations on 
grids over one millions cells. For the largest 
calculations, the SGI Non Uniform Memory 
Architecture (NUMA) high performance computers are 
preferred not only because of performance, but because 
they are true 64 bit computers. (They can address over 
two GB of memory. This is needed for some pre and 
post processing operations.) 

LDV Measurements 

The LDV system used in this paper is described 
in detail by Jeng et al. 11 Only a brief description will be 
given in this paper. The TSI 2D LDV measurement 
system 12 was used to measure the vertical velocity and 
one component of horizontal velocity. To get improved 
signals, forward scatter was used to measure the airflow 
field. The measurement system includes a 3w argon 
laser, an optical fiber driver, a transmitter, a signal 
receiver, a photomultiplier, and a signal processor. The 
optical transmitter and the receiver were mounted on a 
3-D moveable traverse that was controlled by the 
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Figure 6 - LDV Measurement Locations 

computer. A 5 00mm- focus lens transmitter and 300-mm 
focus lens receiver were used in the measurement. The 
measurement involved the mean velocity and turbulent 
characteristics. Though the refraction of the laser beams 
through the walls may lead to the separation of the two 
focus points, the thin walls will decrease this separation. 
By using the coincidence selection on the software, the 
shear stress was obtained. Usually, the data rate of the 
measurements was 600/s and 3000 measurements were 
collected at each data position. The time range is about 5 
seconds per point. 

In the experiment, a TSI six jet aerosol was 
used to generate seeding particles, which have a 
diameter of 3 pm. It almost has no slip relative to the air 
and has good refracting characteristic. Measurements 
were obtained at twelve axial locations downstream of 
the multi-swirler array. Locations of the measurements 
are shown in Figure 6. At each location, due to the 
similarity of the design, only 4/9 of the area was 
measured. The measurement resolution is lxl mm for 
the first five locations and 2x2 mm resolution for the 
other locations. 


Y 



Figure 7 - Anti-Symmetry of the Horizontal 
Velocity Used in the LDV Measurements 


In order to represent the entire flow field, the 
other component of the horizontal velocity, U, was 
obtained by a so called “symmetry method”. The 
method can be explained by Figure 7. In Figure 7, we 
assume that the horizontal velocity is anti-symmetry 
relative to the diagonal line of the area. So we have 
ul=-v2 and vl=-u2. By this method, we obtained the 
information of the other component. Though it is not 
accurate, it gave us enough information to understand 
the horizontal velocity distribution of the flow field. 
The estimated error for the time averaged and RMS 
velocities is 0.25 m/s. 

Data Reduction and Post Processing 

The NCC computations and the LDV 
measurements were both manipulated to give the 
appearance of the full geometry. The NCC calculations 
were periodically mirrored every ninety degrees. This 
resulted in the post processing of over 10 million 
computation cells. For this purpose, a four processor 
SGI Onyx 2 graphical supercomputer was used. The 
Ensight Gold 13 post-processing package was used in 
parallel mode to extract results. The LDV 

measurements were manipulated in a similar fashion. 
The graphical computing requirements for the LDV 
measurements were far less than the NCC 

computations. The Tecplot 14 post-processing package 
was used to manipulate extracted line plot data. 

RESULTS 

The resulting data is compared via XY plots, 
contour and vector plots, and isosurfaces. In addition, 
animations containing additional three-dimensional 
entities like elevated surfaces, stream tubes, and flow 
particles traces were created. These entities will not be 
shown because of paper size restrictions and the fact 
that these entities to do not display well on a two- 
dimensional piece of paper. 

Figure 8 shows isosurfaces for W (or the axial) 
velocity equal to zero. This gives a good overview of 
the two flows fields from both the NCC and LDV data. 
The isosurfaces are similar, but not the same. The NCC 
computations show that the recirculating regions are 
generally shorter that the LDV data. The corner 
swirlers stay centered for the NCC computations, but 
the LDV measurements show that the recirculation 
zones move toward the walls. 

Figure 9 through Figure 14 show combined contour 
plots for velocity magnitude, as well as velocity 
vectors. The contours are on the same scale, but the 
vectors are on a similar, but not the exactly same scale. 
The vectors should only be used to qualitatively 
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compare the NCC and the LDV data. Figure 9 shows 
that at the Z = 3mm plane, both the NCC and LDV data 
produce similar results. The NCC predicts a generally 
stronger recirculation zone away from the jets. For 
Figure 10, Z = 13mm, the NCC shows a more ordered 
flow field, with slightly stronger recirculation zones 
that are more circular. The LDV data shows a more 
distorted flow field, with flow from the corner swirlers 
starting to distort. The LDV data also indicated that the 
jets are still present, while the NCC data smears the jets 
at this location. Figure 11, Z = 23mm, shows that the 
NCC data still has nine strong vortices, while the LDV 
measurements show that most of the strong flow 
features have washed away. Figure 12, Z = 50mm, 
gives similar results, showing that NCC still predicts 
features in the flow, while the LDV data shows that 
most of the flow features have washed away. However, 
these differences between the NCC and LDV data are 
small. 

Figure 13 and Figure 14 are cuts along the X 
plane at 0mm (centerline) and 19mm (between the 
swirlers). Both the NCC and LDV data trends the 
same, by dissipating the jets around 30mm. The NCC 
computations predict stronger flow features that the 
LDV measurements. The LDV measurements are taken 
at 1 and 2mm intervals, which could account for the 
stronger flow feature predicted by the NCC. 

Since the prediction of recirculation zones are 
of primary interest, Figure 15 through Figure 20 show 
contours of the axial (W) velocity. At Z = 3mm, Figure 
15 shows that the NCC and LDV data compare closely. 
The flow is more distorted for the LDV measurements, 
while NCC show circular contours. The NCC produces 
slightly stronger recirculation regions (as given by 
negative axial velocity). Figure 16, Z =13mm, shows 
that the NCC is producing more uniform contours 
compared to the LDV measurements, but the values of 
the contours compare well. Figure 17, Z = 23mm, 
shows that the NCC computations are still producing 
nine strong recirculation regions, while the LDV 
measurements indicate that only four strong 
recirculation zones exist. At Z=50 mm, Figure 18 
shows that flow features wash out for both NCC and 
LDV data. Both show a cross with four vortices, but 
the LDV features are slightly more shifted. Figure 19, 
X = 0mm, show near the same contour lines. Again, 
the scales are different for the LDV measurements and 
the NCC predictions. As previously shown in Figure 8, 
LDV measurements indicate that the center swirler has 
a greater recirculation zone. Figure 20, X = 19mm, the 
LDV data shows only one large recirculation zone, 
while the NCC computations predict that two large 
recirculation zones exist. 




Figure 8 - Isosurface Comparisons between NCC 
Computations and LDV measurements for W 
Velocity = 0 m/s 


(Text continues on page 15) 
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NCC Computations of the LDI 35 Degree Co Rotating Swirler 
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LDV Measurements of the 35 Degree LDI Swirler 
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Figure 9 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the Z = 3mm 
Plane 


Figure 10 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the Z = 13mm 
Plane 


NAS A/TM— 200 1 -2 1 076 1 


6 



NCC Computations of the LDI 35 Degree Co Rotating Swirler 
4 % Pressure Drop, Atmospheric Pressure 

Cubic Non-Linear K-epsilon Turbulence Model, Low Re Wall Treatment 


Evec_Nodal [m/g] 


above 

65.00 

60.00 

55.00 

50.00 

45.00 

40.00 

35.00 

30.00 

25.00 

20.00 

15.00 

10.00 

5.00 

below 







NCC Computations of the LDI 35 Degree Co Rotating Swirler 
4 % Pressure Drop, Atmospheric Pressure 

Cubic Non-Linear K-epsilon Turbulence Model, Low Re Wall Treatment 


Evec_Nodal [m/g] 
above _ 

65.00 ■ 

60.00 ■ 

55.00 

50.00 

45.00 

40.00 

35.00 

30.00 

25.00 

20.00 

15.00 

10.00 U 

5.00 B 
below ■ 



LDV Measurements of the 35 Degree LDI Swirler 
4 % Static Pressure Drop, Atmospheric Pressure 


Velocit y[m/ 5 ] 


above 

65.00 

60.00 

55.00 

50.00 

45.00 

40.00 

35.00 

30.00 

25.00 

20.00 

15.00 

10.00 

5.00 

below 


i 

. 
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Figure 11 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the Z = 23mm 
Plane 


Figure 12 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the Z = 50mm 
Plane 
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Figure 13 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the X = 0mm 
Plane 


Figure 14 - Comparisons between NCC 
Computations and LDV Measurements, Vmag 
Contours and Velocity Vectors at the X = 19mm 
Plane 
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Figure 15 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the Z = 3mm Plane 


Figure 16 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the Z = 13mm Plane 
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Figure 17 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the Z = 23mm Plane 


Figure 18 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the Z = 50mm Plane 
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Figure 19 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the X = 0mm Plane 


Figure 20 - Comparisons between NCC 
Computations and LDV Measurements, Axial (W) 
Velocity Contours at the X = 19mm Plane 
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Figure 21 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at the Z = 3mm 
Plane 


Figure 22 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at the Z = 13mm 
Plane 
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NCC Computations of the LDI 35 Degree Co Rotating Swirler 
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Figure 23 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at the Z = 23mm 
Plane 


Figure 24 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at the Z = 50mm 
Plane 
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Figure 25 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at the X= 0mm 
Plane 


Figure 26 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (K) at the X = 19mm 
Plane 
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At this time, NCC cannot report Reynolds 
stress components. LDV RMS measurements can be 
easily converted to k, the turbulent kinetic energy, by 
Equation (1). The NCC computes k as part of the k - 
epsilon model. 

Equation (1) 



Figure 21 through Figure 26 show contours of turbulent 
kinetic energy. Figure 21, Z = 3mm, shows that NCC 
underpredicts the turbulent kinetic energy by 50 % 
compared to the FDV measurements. The shape of the 
contours is similar. Figure 22, Z = 13mm, the NCC 
calculations underpredict the turbulent kinetic energy, 
and the contours are different. The FDV measurements 
show the turbulent kinetic energy has an even nature, 
which the NCC calculations do not show. Fooking at 
Figure 23, Z = 23mm, the FDV measurements indicate 
that the turbulent kinetic energy does not vary much at 
this plane. NCC computations underpredict the 
turbulent energy and do not show the same contour 
lines. At Z=50mm, Figure 24 show that the turbulent 
kinetic energy, has dissipated for both the NCC 
computations and the FDV measurements. Figure 25 
and Figure 26, X = 0mm and 19mm, show once again 
that NCC calculations underpredict the turbulent energy 
by approximately 50% and generally do not show the 
even nature regarding this turbulent flow. The NCC 
calculations and FDV measurements show that around 
50mm, the turbulent kinetic energy does dissipate to 
nearly zero. 

Contour and Vector plots present good overall 
comparisons but do not quantitatively compare the 
NCC and FDV data. For this reason, line plots are 
extracted from both data sets. Figure 27 shows where 
the line plots were extracted. Figures 28 through 31 
compare velocity magnitude at line plots 1,2, and 3. 
Fine Plot 1 is a line extracted axially (Z) at X = 
25.4mm, Y = -25.4mm (in the center of a corner 
swirler). Fine Plot 2 is a line extracted axially at X = 
12.7mm, Y = -12.7mm (one of four positions between 
the swirl ers). Fine Plot 3 is a XY diagonal line at the Z 
= 23mm plane (at this plane the NCC computations and 
FDV differ). Figure 28 shows a comparison for the 
velocity magnitude for Fine Plots 1-3. At Fine Plot 1, 
the NCC predicts the general trend for the velocity 
magnitude, but at 5 - 10mm and 25 - 35mm, the NCC 
is off by 30 to 50%. For Fine Plot 2, the NCC shows 
the overall trend, but does not show the flow structure 


in the 3 - 10mm region. The NCC generally over 
predicts values by 50% from 10 - 25mm. Past this 
region, the NCC generally agrees within 20% of the 
FDV data. At Fine Plot 3, the NCC computations seem 
to smear out flow structures compared to the FDV 
measurements. The NCC is generally within 20% of 
the FDV data. Figure 29 shows a comparison for the U 
velocity at Fine Plots 1-3. Fine Plot 1 shows that the 
NCC predicts the same flow structure 10mm before the 
FDV measurements. Values are off by more than 50% 
in the 5 - 25mm region, with the NCC agreeing within 
20% after 25m. Fine Plot 2 shows the same general 
trend, disagreeing in the 1 0 - 20mm region. Except for 
this region, the NCC is within 20% of the experimental 
values. At Fine Plot 3, the NCC shows the general 
trend, with values within 20% of the smooth regions. 
In the 22 - 35mm region, a region of a strong gradients, 
the NCC disagrees by more than 50%. Figure 30 shows 
a comparison for the V velocity at Fine Plots 1-3. At 
Fine Plot 1, the trend the NCC predicts is poor. Values 
are off 30 - 50% near the swirler exit and diverge 
downstream. Fine Plot 2 shows a general trend. From 
5 - 40mm , the V velocity is underpredicted 30 - 50%. 
The comparison between the NCC computations and 
FDV measurements gets better past 40mm. For Fine 
Plot 3, a trend is not shown. The NCC values seem to 
be out of phase with the FDV data. Figure 31 shows a 
comparison for the W velocity for Fine Plots 1-3. For 
Fine Plot 1, the agreement is very good. For the entire 
plot, the NCC calculations are within 10 - 20% of the 
FDV measurements. Fine Plot 2 shows good 
agreement within 20% of the FDV data. In the 20 - 
30mm region, the comparison is within 30 -50%. Fine 
Plot 3 shows a general trend, but the NCC does not 
show all of the flow structures. At 15 - 25mm, the 
NCC is off by more than 30 - 50%, and at 35 - 45mm, 
it is off by more than 50%. 
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Line Plot 1 



Line Plot 2 



Line Plot 3 



Line Plot 1 



Line Plot 2 



Line Plot 3 


Figure 28 - Comparisons between NCC 
Computations and LDV Measurements, Velocity 
Magnitude at Line Plots 1, 2, and 3 


Figure 29 - Comparisons between NCC 
Computations and LDV Measurements, U 
Velocity at Line Plots 1, 2, and 3 
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Line Plot 1 



Line Plot 2 



Line Plot 3 



Line Plot 1 



Line Plot 2 



Line Plot 3 


Figure 30 - Comparisons between NCC 
Computations and LDV Measurements, V 
Velocity at Line Plots 1, 2, and 3 


Figure 31 - Comparisons between NCC 
Computations and LDV Measurements, W 
Velocity at Line Plots 1, 2, and 3 
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Line Plot 2 



Figure 32 - Comparisons between NCC 
Computations and LDV Measurements, 
Turbulent Kinetic Energy (k) at Line Plots 1,2, 
and 3 
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Figure 32 shows a comparison for the turbulent kinetic 
energy (k) for Line Plots 1-3. At Line Plot 1, the 
agreement is good, with the value of k within 20% of 
the experimental data. Line Plot 2 shows a general 
trend, but the value is off by more than 50% in the 5 - 
25mm region. For Line Plot 3, a general trend is 
shown, with data off by more than 50% in the 0 - 
20mm region. Past 20mm, the NCC agrees with the 
LDV data to within 20 - 30% (approaching the corner). 

DISCUSSION 

The NCC computations tend to compare well 
with the LDV measurements in some regions and 
compare poorly in other regions of the LDI swirler. 
The NCC computations generally show the same trends 
as the experimental data. However, the LDV data 
cannot be considered the absolute reference at this time. 
Because of the symmetry method, the U and V 
velocities are essentially the same (they are opposite). 
The U velocity is the measured velocity component. 
The NCC computations diverge for the V velocity 
downstream of the injector. This trend is not consistent 
with any other variable. Because of this, the V velocity 
should be measured. The W Velocity compares well 
for most of the one, two, and three-dimensional 
representations. While the turbulent kinetic energy is 
not the same for all of the contour plots, the line plots 
indicate a trend. The turbulent kinetic energy is also 
affected by the measurement. (See Equation (1) 
for the definition of k.) The NCC can predict the 
general trends for the vast majority of the flow field, 
and compares well in specific regions. Because of this, 
the NCC can be used as a design tool. Nevertheless, 
improvements should be made to the NCC. 

The line plots indicate that the NCC poorly 
predicts flow quantities in regions of strong gradients. 
To improve this, two steps should be taken. The first 
step is to scale the numerical dissipation with matrix 
dissipation instead of scalar numerical dissipation. The 
methods suggested by Swanson and Turkel for non - 
reacting flows 15 , and Gerlinger et al for reacting flows 
should be used. 16 Results indicate that matrix 
dissipation greatly improves the accuracy of central 
difference finite volume CFD codes. The second step is 
to add Adaptive Mesh Refinement (AMR) to the NCC. 
This is needed because it is nearly impossible to 
correctly setup a computational mesh for very complex 
flows. The Pyramid 3D package is an excellent 
candidate to do the AMR in NCC. 17 

Turbulence modeling also needs to be 
improved. A “unified” turbulence model should be 
implemented in the NCC. 18 This model would combine 
a wall distance free, cubic non-linear, k-epsilon, low 
Reynolds turbulence model with the current dynamic 
wall function. This would give the benefits of a low 


Reynolds number model in regions of high gradients 
without the sensitivity of grid point placement. In 
addition, the Reynolds stresses should be stored during 
NCC computations. 

CONCLUSIONS 

A Comparison between NCC computations 
and LDV measurements was made. This comparison 
showed that the NCC was able to predict the general 
trends of the flow field. Because of this study, the NCC 
can be used as a design tool for LDI injectors. 
Suggestions were given on how to further improve the 
NCC for turbulent, three-dimensional, recirculating 
flows. 
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